SYNCHRONIZATION OF SPECTRAL COMPONENTS 
AND ITS REGULARITIES IN CHAOTIC DYNAMICAL SYSTEMS * 
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The chaotic synchronization regime in coupled dynamical systems is considered. It has been 
shown, that the onset of synchronous regime is based on the appearance of the phase relation 
between interacting chaotic oscillators frequency components of Fourier spectra. The criterion of 
synchronization of spectral components as well as the measure of synchronization have been dis- 
cussed. The universal power law has been described. The main results are illustrated by coupled 
Rossler systems, Van-der-Pol and Van-der-Pol-Duffing oscillators. 
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Introduction 

Chaotic synchronization is one of the fundamental phe- 
nomena actively studied last time yj, having impor- 
tant both theoretical and applied significance (e.g., for 
the information transmission by means of deterministic 
chaotic signals 0,0], in biological Q and physiological Q 
tasks, etc.). Several different types of chaotic synchro- 
nization of coupled oscillators, i.e. generalized synchro- 
nization , phase synchronization [3, lag synchroniza- 
tion and complete synchronization Q are traditionally 
distinguished. There are also attempts to find a unifying 
framework for chaotic synchronization of coupled dynam- 
ical systems 0, CU Ull. 

In our works [T^ . Il3| it was shown that phase, gener- 
alized, lag and complete synchronization are closely con- 
nected with each other and, as a matter of fact, they 
are different manifestations of one type of synchronous 
oscillations behavior of coupled chaotic oscillators called 
time scale synchronization. Synchronous regime charac- 
ter (phase, generalized, lag or complete synchronization) 
is defined by the presence of synchronous time scales 
s, introduced by means of continuous wavelet trans- 
form 0, ITU ITfj| with Molrlet mother wavelet function. 
Each of time scales can be characterized by the phase 
4> s (t) — arg W(s, t), where W(s, t) is the complex wavelet 
surface. In this case, the phenomenon of chaotic syn- 
chronization of coupled systems is manifested by a syn- 
chronous behavior of the phases of coupled chaotic os- 
cillators 4>si,2{t) observed on certain synchronized time 
scales range s m < s < s&, for time scales s from which 
the phase locking condition 



\4>si(t) - <fi S 2(t)\ < const 



(1) 



is satisfied, and the part of the wavelet spectrum ene rgy 
being fallen on this range does not equal zero (see [T3 . Tl7| 
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for detail). The range of synchronized time scales [sm 5 Sb\ 
expands when the coupling parameter between systems 
increases. If the coupling type between oscillators is 
defined in such a way that the lag synchronization ap- 
pearance is possible then all time scales become synchro- 
nized with further coupling parameter increasing, while 
the coinciding states of interacting oscillators are shifted 
in time relative to each other: xi(t — r) ~ X2(t). Fur- 
ther coupling parameter increasing leads to a decrease of 
the time shift r. The oscillators tend to the regime of 
complete synchronization, Xi(t) ~ ^(t) and the phase 
difference (cf) s i(t) — <f> S 2(t)) tends to be zero on all time 
scales. 

The time scale s introduced into consideration by 
means of continuous wavelet transform can be consid- 
ered as a quantity which is inversely proportional to the 
frequency / defined with the help of Fourier transforma- 
tion. For the Morlet mother wavelet function |lf| with 
parameter f2 = 2ir the relationship between the frequency 
/ and the time scale is quite simple: s = 1//. Therefore, 
the time scale synchronization should also manifest in 
the appearance of the phase relation between frequency 
components / of corresponding Fourier spectra S(f) of 
interacting oscillators. 

In this paper we consider the synchronization of spec- 
tral components of Fourier spectra of coupled oscilla- 
tors. We discuss the mechanism of chaotic synchroniza- 
tion regime manifestation in coupled dynamical systems 
based on the appearance of the phase relation between 
frequency components of Fourier spectra of interacting 
chaotic oscillators (see also 18]). One can also consider 
the obtained results as a criterion of existence (or other- 
wise, an impossibility of existence) of lag synchronization 
regime in coupled dynamical systems. 

The structure of this paper is the following. In Sec- 
tion [I] we discuss the synchronization of spectral compo- 
nents of Fourier spectra. We illustrate our approach with 
the help of two coupled Rossler systems in Sec. |H] The 
quantitative measure of synchronization is described in 
Section IlIII The universal power law taking place in the 
presence of time scale synchronization regime is discussed 
in Sections IIVI and Ivl The final conclusion is presented 
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in Section IVT1 

I. SYNCHRONIZATION OF SPECTRAL 
COMPONENTS OF FOURIER SPECTRA 

It should be noted that the continuous wavelet trans- 
form is characterized by the frequency resolution lower 
than the Fourier transformation (see [Tit EJ). The 
continuous wavelet transform appears as smoothing of 
the Fourier spectrum, whereby the dynamics on a time 
scale s is determined not only by the spectral compo- 
nent / = 1/s of the Fourier spectrum. This dynamics is 
also influenced by the neighboring components as well, 
the degree of this influence depends both on their po- 
sitions in the Fourier spectrum and on their intensities. 
Thus, the fact that coupled chaotic oscillators exhibit 
synchronization on a time scale s of the wavelet spectrum 
by no means implies that the corresponding components 
/ = 1/s of the Fourier spectrum of these systems are also 
synchronized. 

Let Xi (t) and x 2 (t) be the time series generated by the 
first and the second coupled chaotic oscillators, respec- 
tively. The corresponding Fourier spectra are determined 
by the relations 

+00 

Si,a(/)= / zi, 2 (i)e- ,;27r/ * dt. (2) 

—00 

Accordingly, each spectral component / of the Fourier 
spectrum S(f) can be characterized by an instantaneous 
phase 4>f(t) = 0/o + 2irft, where 0/o = axgS(f). How- 
ever, since the phase 0/(i) corresponding to the fre- 
quency / of the Fourier spectrum S(f) increases with the 
time linearly, the phase difference of the interacting os- 
cillators at this frequency 0/i(i) — 0/2(0 = 0/oi ~ 4>f 02 
is always bounded and, hence, the traditional condition 
of phase entrainment (used for detection of the phase 
synchronization regime) 

|0i(t) -0 2 (t)| < const, (3) 

is useless. Apparently, a different criterion should be used 
to detect the synchronization of coupled oscillators at a 
given frequency /. 

In the regime of lag synchronization, the behavior of 
coupled oscillators is synchronized on all time scales s of 
the wavelet transform (see |l2j). Therefore, one can ex- 
pect that all frequency components of the Fourier spectra 
of the systems under consideration should be synchro- 
nized as well. In this case, x\{t — r) ~ x 2 (t) and, hence, 
taking into account J5J) one has to obtain 

S 2 (f) ~ S 1 (f)e l2 ^f. (4) 



Thus, in the case of coupled chaotic oscillators occurring 
in the regime of lag synchronization their instantaneous 
phases corresponding to the spectral component / of the 
Fourier spectra Si^(f) will be related to each other as 
4>f 2 (t) — <t>fi(t) + 2m~f and, hence, the phase difference 
0/2 (t) — 4>fi(t) of coupled oscillators on the frequency / 
must obey the relation 

A0 / = /1 (i)-0 /2 (t) = 

= 0/oi - 0/02 = 27TT/. 

Accordingly, the points corresponding to the phase differ- 
ence Aipf of the spectral components of chaotic oscilla- 
tors in the regime of lag synchronization on the (/, A0/) 
plane must fit a straight line with slope k — 2ttt. In the 
case of complete synchronization of two coupled identi- 
cal oscillators the slope of this line k is equal to zero (see 
also [13). 

The destroying of the lag synchronization regime (e.g., 
as a result of decrease of the coupling strength between 
oscillators) and transition to the regime of phase syn- 
chronization (in the case when the instantaneous phase 
of chaotic signal can be introduced correctly [2(j ) results 
in the loss of synchronism for a part of time scales s of the 
wavelet spectra ^2j • Accordingly, one can expect that a 
part of spectral components of the Fourier spectra in the 
phase synchronization regime will also lose synchronism 
and the points on the (/, A0/) plane will deviate from the 
straight line mentioned above [3{|. It is reasonable 
to assume that synchronism will be lost primarily for the 
spectral components / characterized by a small fraction 
of energy in the Fourier spectra Si^(f), while the com- 
ponents corresponding to a greater energy fraction will 
remain synchronized and the corresponding points on the 
(/, A0/) plane will be located at the straight line as be- 
fore. As the lag synchronization regime does not occur in 
the system any more, the time shift r can be determined 
by the delay of the most energetic frequency component 
f m in the Fourier spectra t = (0/ m 2 — 0/ m i)/(27r/ m )- 

As the coupling parameter decreases further, an in- 
creasing part of the spectral components will deviate 
from synchronism. However, as long as the most "en- 
ergetic" components remain synchronized, the coupled 
systems will exhibit the regime of time scale synchro- 
nization. Obviously, that for the synchronized spectral 
component the phase difference At/?/ is located after the 
transient finished independently on initial conditions. 

To describe the synchronization of spectral compo- 
nents, let us introduce a quantitative characteristic of 
a number of spectral components of the Fourier spectra 
<Si,2(/) occurring in the regime of synchronism, 
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+ 00 



CL = 



/ ff(|5i(/)| 2 - L)H(\S 2 (f)\ 2 - L) (A<f> f - 2vTf) z df 
_o 

+00 

I H(\Si(f)\2-L)H(\S 2 (f)\2-L)df 




(6) 



r 



where H(£) is the Heaviside function, L is the threshold 
power level (in dB) above which the spectral components 
are taken into account, and r is determined by the time 
shift of the most energetic frequency component (f m ) 
in the Fourier spectra, r = (4>f m 2 — <^/ m i)/(27r/ TO ). The 
quantity <tl tends to be zero in the regimes of complete 
and lag synchronization. After the lag synchronization 
regime destroying caused by the decrease of the coupling 
strength the value of ctl increases with the number of 
desynchronized spectral component of the Fourier spec- 
tra S\fi{f) of coupled oscillators. 

Real data are usually represented by discrete time se- 
ries of a finite length. In such cases, the continuous 
Fourier transform (0) has to be replaced by its discrete 
analog (as it was done in \Vt^), and the integral by 
the sum 



1 N 



(7) 



taken over all spectral components of the Fourier spectra 
<Si,2(/) with the power above L. In calculating ctl, it is 
expedient to perform averaging over a set of time series 
£1,2 (i)- The phase shift A<pf can be calculated either 
in the wa y as it was done in or by means of cross 
spectrum [21j . 



II. TWO MUTUALLY COUPLED ROSSLER 
SYSTEMS SYNCHRONIZATION 

In order to illustrate the approach proposed above, let 
us consider two coupled Rossler systems 

Xl,2 = -U\.2V\.2 - 21,2 + £(^2,1 ~ Zl,2), 

2/1,2 = ^1,2^1,2 + ayi,2 + £(2/2,1 - 2/1,2), (8) 
ii,2 = P + 21,2(2:1,2 - c), 

where e is the coupling parameter, io\ = 0.98, L02 = 1.03. 
By analogy with the case studied in |22|, the values 
of control parameters have been selected as follows: 
a = 0.22, p = 0.1 and c = 8.5. It is known [H that 
two coupled Rossler systems with e = 0.05 occur in the 
regime of phase synchronization, while for e = 0.15 the 
same systems exhibit lag synchronization. 

Figure l,a shows a plot of the value ox versus cou- 
pling parameter e. One can see that ox tends to be zero 
when the coupling parameter e increases, which is the 



evidence of the transition from phase to lag synchroniza- 
tion. Figures 1,6—1,/ illustrate the increase in the num- 
ber of synchronized spectral components of the Fourier 
spectra S\ji{f) of two coupled systems with coupling 
strength e increasing. Indeed, Fig. 1,6 corresponds to the 
asynchronous dynamics of coupled oscillators (e = 0.02). 
There are no synchronous spectral components for such 
coupling strength and dots are scattered randomly over 
the (fjAtfif) plane. The weak phase synchronization 
(e = 0.05) after the regime occurence is shown in the 
Fig. l,c. There are a few synchronized spectral com- 
ponents the phase shift Acff of which satisfies the con- 
dition Almost all spectral components are non- 
synchronized, therefore the points corresponding to the 
phase differences Aipj are spread over (/, A<fif) plane and 
the value of <jl is rather large. 

Figures l,d and 1, e correspond to the well pro- 
nounced phase synchronization (e = 0.08 and 0.1, respec- 
tively). Fig. 1,/ shows the state of lag synchronization 
(e = 0.15), when all spectral components / of Fourier 
spectra are synchronized. Accordingly, in this case all 
points on the (f,A<j>f) plane are at the line © with 
slope k = 27tt. With coupling strength e increasing, the 
value of ox decreases monotonically that verifies the as- 
sumption that when two coupled chaotic systems undergo 
transition from asynchronous dynamics to lag synchro- 
nization, more and more spectral components become 
synchronized. When all spectral components / are syn- 
chronized, the phase shift A(ff for them is 27rr/, there- 
fore the points on the (f,A(f>f) plane lie on the straight 
line 10 and the value of ol is equal to zero. 

Another important question is which spectral compo- 
nents of the Fourier spectra of interacting chaotic oscilla- 
tors are synchronized first and which do it last. Figure 2, a 
shows a plot of the ol value for the coupling strength 
e = 0.05 (corresponding to the weak phase synchroniza- 
tion) versus power L at which the spectral components 
fj of the Fourier spectra Si,2(/) are taken into account 
in equation (JJJ). One can see that the "truncation" of 
the spectral components with small energy leads to a 
decrease of the ol value. Figures l,6-l,e illustrate the 
distribution of the phase difference Aipf of the spectral 
components / with the power exceeding the preset level 
L. The data in Fig. 2 show that the most "energetic" 
spectral components are first synchronized upon the on- 
set of time scale synchronization. On the contrary, the 
components with low energies are the first to go out from 
synchronism upon destroying of the lag synchronization 
regime. 
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FIG. 1: (a) The value <jl versus coupling parameter e and 
(&)-(/) the phase difference Aipf of various spectral compo- 
nents / of the Fourier spectra Si,2(/) of two coupled Rossler 
systems for different values of coupling strength e. (b) The 
asynchronous dynamics for coupling parameter e = 0.02; (c) 
The chaotic synchronization regime e = 0.05; (d), e = 0.08; 
(e) e — 0.1 and (/) e = 0.15. The plots are constructed for the 
time series asi,a(i) with a length of 2000 dimensionless time 
units at a discretization step of h = 0.2 at a power level of 
L = —40 dB of Fourier spectra Sipif) 



III. A CRITERION AND MEASURE OF 
SYNCHRONIZATION 

Let us briefly discuss a criterion of spectral components 
synchronization. Obviously, the relation (J5J is quite con- 
venient as a criterion of synchronism in case of lag syn- 
chronization destroying in the way considered above. If 
the type of coupling between systems has been defined 
in such a manner that the lag synchronization regime 
can not appear, the relation JSJ can not be the criterion 
of spectral components synchronization. So, as a general 
criterion of synchronism of identical spectral components 
/ of coupled systems we have to select the different con- 
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FIG. 2: (a) The value ox versus power L at which the spectral 
components fj of the Fourier spectra Si,2(/) are taken into 
account in equation Q. (b)-(e) The phase difference A</?/ of 
various spectral components / of the Fourier spectra £>i,2(/) 
of two coupled Rossler systems for various power levels L = 
-40 dB (6), L = -30 dB (c), L = -20 dB (d) and L = -10 
dB (e) for coupling strength e = 0.05 



dition rather than (0 . As such a criterion we have chosen 
the establishment of the phase shift 



Aipf = 0/oi - 0/02 = const 



(9) 



which must not depend on initial conditions. To illustrate 
it let us consider the distribution of phase difference Atp t 
obtained from the series of 10 3 experiments for Rossler 
systems (JSJ. Fig. 3, a corresponds to the asynchronous dy- 
namics of coupled oscillators when the coupling param- 
eter e — 0.02 is below the threshold of chaotic synchro- 
nization appearance (see also Fig. l,a). One can see, that 
the phase difference Aipf for the spectral components / 
of Fourier spectra Si^t) in this case is distributed ran- 
domly over all interval from — n to n. It means that the 
phase shift between spectral components / is different for 
each experiment (i.e., for different initial conditions) and, 
therefore, there is no synchronism whereas the considered 
frequency / is the same for both spectra Si^(f)- The 
similar distributions are observed for all spectral com- 
ponents / in the case of asynchronous regime (see also 
Fig. l,o and Fig. 3,&), though one can distinguish the 
maximum in the distribution on the frequency / close 
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FIG. 3: Distribution N(Acj)f) of phase difference Ay;^ ob- 
tained from series of 10 3 experiments for Rossler systems (|HJ . 
(a) The asynchronous dynamics takes place (e = 0.02), the 
distribution of phase shift Aipf has been obtained for the 
spectral components / ~ 0.0711; (6) e = 0.02, / ~ 0.1764; 
(c) the distribution of phase shift for synchronous spectral 
component / ~ 0.1764 (e = 0.08) and (d) the analogous dis- 
tribution for asynchronous spectral component / ~ 0.0711 for 
the same coupling parameter e = 0.08. Compare with Fig. 1 



to the main frequency of Fourier spectrum S(f) as an 
prerequisite of synchronization beginning. 

When the systems demonstrate synchronous behavior 
the distributions N(Acj)f) of phase shift A<j>f are quite 
different. In this case one can distinguish both synchro- 
nized and non-synchronized spectral components charac- 
terized by distributions of phase shift of different types. 
In Fig. 3,c the distribution of A<j>f for synchronous spec- 
tral component A(ff is shown. One can see that it looks 
like 5-function that means the phase shift Acfit is always 
the same after transient finished. Obviously, this phase 
shift A<fif does not depend on initial conditions. 

For the non-synchronized spectral components the dis- 
tributions N(A(j)f) are different (see Fig. 3,rf). Evidently, 
in this case the phase shift A<f)f is varied from experi- 
ment to experiment. At the same time, the tendency to 
the synchronization of these spectral components can be 
observed. The distribution N(A<fif) looks like Gaussian. 
With the coupling parameter increasing the dispersion of 
it decreases and the spectral components / of considered 
Fourier spectra Si^(f) tend to be synchronized. The 
same effect can be observed in Fig 1,6-/. With coupling- 
parameter e increasing, the points on the (/, A<fif) plane 
tend to fit a straight line with slope k = 2t:t and their 
scattering decreases. 

So, the general criterion of synchronism of identical 




FIG. 4: The dependence of the synchronization measure p for 
the first Rossler system © on coupling parameter e 



spectral components / of coupled systems is the estab- 
lishment of the phase shift © after transient finished. 
It is important to note that the case of classical synchro- 
nization of periodical oscillations also obeys to considered 
criterion © (see, e.g., p3|l- 

Let us consider now the quantitative characteristic of 
synchronization. In [12] the measure of synchroniza- 
tion based on the normalized energy of synchronous time 
scales has been introduced. The analogous quantity p 
may be defined for Fourier spectra S(f) as 



Pi,a = p J \SiAf)\ 2 df 



(10) 



where F s is the set of synchronized spectral components 
and 



+ OG 



P 



\SiAf)\ 2 df 



(11) 



is the full energy of chaotic oscillations. In fact, the value 
of p is the part of the full system energy corresponding 
to synchronized Fourier components. This measure p is 
for the nonsynchronized oscillations and 1 for the case 
of complete and lag synchronization regimes as well as 
the quantity introduced in 12]. When the systems un- 
dergo transition from asynchronous behavior to the lag 
synchronization regime the measure of synchronism takes 
the value between and 1 that corresponds to the case 
when there are both synchronized and nonsynchronized 
spectral components in Fourier spectra Si^if)- 

For the real data represented by discrete time series 
of a finite length one have to use the discrete analog of 
Fourier transform while the integrals in the relation llOjl 
should be replaced by the sums 



Pl,2 



A/, 



(12) 
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where 



P = ^|5 1)2 (/ j )| 2 A/. 



(13) 



one can obtain the equations for amplitudes r± t2 and 
phases tpxp of coupled oscillators as follows: 



While the sum in the equation l|12|) is being calculated 
only synchronized spectral components fj should be 
taken into account. 

Figure 4 presents the dependence of the synchroniza- 
tion measure p for the first Rossler oscillator of system JHJ 
on the coupling parameter e. It is clear that the part of 
the energy corresponding to the synchronized spectral 
components grows with the coupling strength increasing. 



IV. SPECTRAL COMPONENTS BEHAVIOR AT 
THE PRESENCE OF SYNCHRONIZATION 

Let us now consider how the closed frequency compo- 
nents of two coupled oscillators behave with increase of 
the coupling strength e. As a model of such situation let 
us select two mutually coupled Van-der-Pol oscillators 



ha = r (A - |?-i,2| 2 ) r ia ± e -^- sm(y>i, 2 - <^2,i), 



<Pl,2 = ' n T -Z — COS(^i^2 - V?2,l)- 



(18) 



2^ 



The oscillations of two Van-der-Pol generators l|14|) are 
synchronized when conditions 



»%2 = 0, lpi 2 = 



(19) 



are satisfied. Assuming the phase difference of oscilla- 
tions Ap = if2 — pi is small enough and taking into 
account only components of first infinitesimal order over 
Ap, one can obtain the relation for the phase shift 



X\ tt 2 + A 2 ± 2JnA(e + QA) 
A^i, 2 = ^-4tt (20) 



2e + 4f7A 



051,2 



(\ - x 2 ia ) x lt2 + tt 2 12 x lt2 = ±e(x 2 ,i -2:1,2), (14) and frequency 



where Q\ j2 = ± A are slightly mismatched natural 
cyclic frequencies, 2:1,2 are the variables, describing the 
behavior of the first and the second self-sustained oscil- 
lators, respectively. The parameter e characterizes the 
coupling strength between oscillators. Nonlinearity pa- 
rameter A = 0.1 has been chosen small enough in order 
to make oscillations of self-sustained generators close to 
the single frequency ones. 

Asymmetrical type of coupling in system l|14|) ensures 
the appearance of the synchronous regime which is sim- 
ilar to the lag synchronization in chaotic systems. For 
such type of coupling the oscillations in the synchronous 
regime are characterized by one frequency w = 2irf while 
small phase shift Atpf between time series Xi t2 (t), de- 
creasing when coupling strength increases, takes place. 

Using the method of complex amplitudes, the solution 
of equation (|14fl can be found in the form 



wi,2 = V^ 2 + A 2 ± 2y/QA(e + QA), 



(21) 



which correspond to the stable and nonstable fixed points 
of the system (11811 . From relations (|2()|) and (|21|l one can 
see that the phase difference Aip of coupled generators is 
directly proportional to the frequency of oscillations uj 
and inversely proportional to the coupling parameter e 
for the small values of detuning parameter A 



Xuj 

^~ 2e" 



(22) 



So, in the synchronous regime the phase shift Ap for 
synchronized frequencies obeys the relation 



Ap 



(23) 



X1.2 = Axoe 1 ' 



A*e 



It is important to note, that the time delay between syn- 
chronized spectral components 



+ A\ 2 e- iwt = 0, 



(15) 



A^ 

U! 



(24) 



where "*" means complex conjugation, uj is the cyclic fre- 
quency, at which oscillations in system (|14f> are realized. 
One can reduce equations (|14|) and i|15|) to the form 

A l , 2 = \{\-\A\ 2 )A+ 

l h ^ 2 ~ ^ Aia T £{A2A ~ Ai ' 2) } (16) 

by means of the averaging over the fast changing vari- 
ables. 

Choosing the complex amplitudes in the form of 



^1,2 = r 1<2 e 



Vi, 2 



(17) 



does not depend upon the frequency, and therefore, the 
time delays for all frequencies / are equal to each other. 
Accordingly, the phase shift Apf for the frequency / 
obeys the relation |JSJ| which is the necessary conditions 
for lag synchronization appearance. Evidently, if the type 
of coupling between oscillators is selected in such manner 
that the phase shift Atpf of synchronized spectral com- 
ponents satisfies the condition l|23|l . the appearance of 
the lag synchronization regime is possible for the large 
enough values of the coupling strength. Otherwise, if the 
established phase shift does not satisfy the condition 12311 , 
the realization of the lag synchronization regime in the 
system is not possible for such kind of coupling. So, the 
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tions 



FIG. 5: The dependence of time shift r between base spec- 
tral components (■ symbols) on coupling parameter e for two 
coupled Rossler systems JHJ . The straight line corresponds to 
the power law r ~ e -1 . The value of coupling parameter 
Ei ~ 0.14 corresponding to the appearance of the lag synchro- 
nization regime is shown by the arrow 



relation (|23[1 can be considered as the criterion of the 
possibility of the existence (or, otherwise, impossibility 
of existence) of lag synchronization regime in coupled 
dynamical systems. 

The regularity (|24[1 takes place for the large number 
of dynamical systems and, probably, is universal. Let us 
consider the manifestations of this regularity for several 
examples of coupled chaotic dynamical systems. 

As the first example we consider the coupled Rossler 
systems (JSJ described above. Obviously, one has to con- 
sider the phase shift A(ff (or time shift r) of synchro- 
nized spectral components to verify the relation l|24|) . As 
it has been shown above, spectral components charac- 
terized by the large value of energy become synchronized 
first when coupling strength increases. So, the main spec- 
tral components f m of Fourier spectra of coupled systems 
are synchronized in the most lengthy range of the coupled 
parameter value. Therefore, it is appropriate to consider 
the time shift r of main spectral components for coupling 
strength values e > 0.05. 

In Fig. 5 the dependence of time lag r between Fourier- 
spectra base frequency components of interacting chaotic 
oscillators on coupling parameter e is shown. Base fre- 
quency u) m = 27r/ m of spectrum is close to u) = 1 and 
slightly changes with coupling parameter increasing. In 
Fig. 5 one can see that after entrainment of Fourier- 
spectrum base spectral components of interacting oscilla- 
tors (that corresponds to establishment of time scale syn- 
chronization regime, see also ^3|) the time lag r, which 
is between them, obeys the universal power law l|24|). 

As the second example we consider the chaotic syn- 
chronization of two unidirectionally coupled Van-der- 
Pol-Duffing oscillators 0, 0, • The drive generator 
is described by system of dimensionless differential equa- 



±i = -v\ \x\ - axi - yi] , 
yi = X\ - j/i - zi, 
z\ = Pyi, 



(25) 



while the behavior of the response generator is defined 
by the system 

±2 = —vi \x\ - ax 2 - 2/2] + V2£(xi - x%), 

2/2 = %2 - V2 - z 2 , (26) 

Z2 = A/2, 

where £1,2, 2/1,2 and z%,2 are dynamical variables, char- 
acterizing states of the drive and response generators, 
respectively. Values of control parameters are chosen as 
following: a = 0.35, (3 = 300, v x = 100, v 2 = 125, the 
difference of parameters 1^1 and z/2 provides the slight 
nonidentity of considered generators. 

In Fig. 6 the dependence of time lag r between time 
realizations of coupled oscillators on coupling parameter 
value e is shown. In this range of coupling parameter 
values the lag synchronization regime is realized. Obvi- 
ously, the time lag r also obeys the power law r ~ e n with 
exponent n = — 1, that corresponds to the relation 124|l . 



V. UNSTABLE PERIODIC ORBITS 

It is important to note another manifestation of power 
law (|24p. It is well known that the unstable periodic 
orbits (UPOs) embedded into chaotic attract ors p lay the 
important role in the system dynamics p l l27 l l28 | includ- 
ing the cases of chaotic synchronization regimes p9l I30L 
l3ll | . The synchronization of two coupled chaotic systems 




FIG. 6: The dependence of time shift r between time series 
xi(t) and X2(t) (■ symbols) on coupling parameter e for two 
unidirectionally coupled chaotic oscillators 125H and 1261 . The 
straight line corresponds to the power law r ~ e _1 



8 




2 4 6 8 10 t 

b 

FIG. 7: (a) The unstable periodic orbit of the length m = 2 
embedded in the chaotic attractor of the first system, (b) 
Time series xi(t) and X2{t) corresponding to the "in-phase" 
unstable saddle orbits of length m — 2 in the first (solid line) 
and the second (dashed line) Rossler systems, respectively. 
The coupling parameter is chosen as e — 0.07. Time shift r 
is denoted by means of the arrow 



in terms of unstable periodic orbits has been discussed in 
detail in ^ij^. It has been shown that UPOs are also syn- 
chronized with each other when the chaotic synchroniza- 
tion in the coupled systems realizes [3^|. Let us consider 
the synchronized saddle orbits m : n (m = n = 1, 2, . . . ), 
where m and n are the length of unstable periodic orbits 
of the first and the second Rossler systems (jHJ), respec- 
tively. It was shown that such synchronized orbits may 
be both "in-phase" and "out-of-phase" , but only "in- 
phase" orbits exist in all range of coupling parameter 
values starting from point of the synchronization begin- 
ning (see for detail). It is known that the time shift 
between synchronized "in-phase" orbits decreases with 
coupling strength increasing. As the UPOs have an influ- 
ence on the system dynamics (and on the Fourier spectra 
of the considered systems, too), it seems to be interest- 
ing to examine whether the time shift r between UPOs 
obeys the power law i|24|). 

To calculate the synchronized saddle orbits we have 
used the SD-method [33], in the same way as it had 
been done in . The UPO embedded in the chaotic 
attractor of the first Rossler system and the time se- 
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FIG. 8: The dependence of time shift r between time series 
xi(t) and X2(t) corresponding to the synchronized saddle or- 
bits of the first and the second Rossler systems on coupling 
parameter e. The symbols □ correspond to unstable orbits 
with length m = 1, the symbols O demonstrate the time shift 
r for the orbits with length m — 2 and the symbols ■ show 
this dependence for orbits with length m = 3. The straight 
line corresponds to the power law r ~ e -1 

ries Xi t 2(t) corresponding to the "in-phase" synchro- 
nized UPOs realized in system (JSJ for coupling strength 
e = 0.07 is shown in Fig. 7. One can see the presence of 
the time shift r between these time series which can be 
easily calculated. 

The calculated time shift r between such "in-phase" 
synchronized saddle orbits appears to obey the power law 
t ~ e -1 as well as the spectral components of Fourier 
spectra do (see Fig. 8). We have examined this relation 
for "in-phase" UPOs with length m = 1 -j- 6 and found 
the time shift dependence on coupling strength agrees 
with power law (|24|l well, but the date is shown in Fig. 8 
only for UPOs with length m = 1 -j- 3 for the clearness 
and the simplicity. So, the power law l|24|) seems to be 
universal and is manifested in different ways. 

VI. CONCLUSION 

In conclusion, we have considered the chaotic synchro- 
nization of coupled oscillators by means of Fourier spec- 
tra; several regularities have been observed. 

The chaotic synchronization of coupled oscillators is 
manifested in the following way. Starting from the cer- 
tain coupling parameter value the synchronization of the 
main spectral components of Fourier spectra of interact- 
ing chaotic oscillators takes place. Therefore, for these 
spectral components / the condition J§J) is satisfied. In 
this case one can detect the presence of the time scale 
synchronization regime (see [l2l Il7|). If for the consid- 
ered systems one can introduce correctly the instanta- 
neous phase of chaotic signal 0, [2(j , one will also detect 
easily the phase synchronization by means of traditional 
approach (see, e.g., 0). 

With further coupling parameter increasing, more 
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spectral components become synchronized. If the cou- 
pling between interacting systems is selected in a such 
way that the lag synchronization regime can be realized 
then the time shift between synchronized spectral com- 
ponents obeys the power law l|24[l . The spectral compo- 
nents characterized by the large value of energy become 
synchronized first. Accordingly, the part of energy being 
fallen on the synchronized spectral components increases 
from (asynchronous dynamics) to 1 (the lag synchro- 
nization regime). Synchronization of all frequency com- 
ponents corresponds to the appearance of the lag syn- 
chronization regime. With further coupling strength in- 
crease, the time lag r obeying the relation (|24|l tends to 
be zero, and related oscillations tend to demonstrate the 
complete synchronization regime. In this case the time 
shift r between synchronized components does not de- 
pend on frequency / of considered components (it is the 
same for all synchronized components) and obeys the 
power law Q24J) with the exponent n = —1. The time 
shift between synchronized "in-phase" UPOs embedded 
in chaotic attractors also obeys the same power law. 

So, in the present paper the mechanism of chaotic 
synchronization regime appearance in coupled dynami- 
cal systems, based on the arising of the phase relation 



between frequency components of Fourier spectra of in- 
teracting chaotic oscillators has been discussed. The ob- 
tained results concerning the power law (|24J) may be also 
considered as the criterion of possible existence (or, oth- 
erwise, impossibility of existence) of lag synchronization 
regime in coupled dynamical systems (i.e., the lag syn- 
chronization regime can not be observed in the coupled 
chaotic oscillators system unless the time shift between 
synchronized components obeys the power law H24|) L 
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